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2nd-order Fermi acceleration as the origin of the Fermi bubbles 
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Abstract: Gamma-ray data from Fermi-LAT show a bi-lobular structure extending up to 50 degrees above and below 
the Galactic centre, coincident with a possibly related structure in the ROSAT X-ray map which presumably originated 
in some energy release close to the centre a few million years ago. It has been argued that the gamma-rays arise due to 
inverse Compton scattering of relativistic electrons accelerated at plasma shocks present in the bubbles. We explore the 
alternative possibility that the relativistic electrons undergo stochastic 2nd-order Fermi acceleration in the entire volume 
of the bubbles by plasma wave turbulence. This turbulence is generated behind the outer shock and propagates into the 
bubble volume, leading to a non-trivial spatial variation of the electron spectral index. Rather than a constant volume 
emissivity as predicted in other models we find an almost constant surface brightness in gamma-rays and also reproduce 
the observed sharp edges of the bubbles. We comment on possible cross-checks in other channels. 
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1 Introduction 

Recently, data from the Fermi-LAT have revealed Q] [Z) 
(see also 0) the presence of two huge bi-lobular structures 
in gamma-rays, the so-called "Fermi bubbles", extending 
up to 50° above and below the galactic plane. The overall 
spectrum of the bubbles is oc E~ 2 , i.e. much harder than 
the 7T°, inverse Compton (IC) and bremsstrahlung fore- 
grounds from galactic cosmic rays in the disk, and extends 
from a spectral shoulder at about a GeV up to a cut-off/roll- 
over at a few hundreds of GeV. The bubbles have an almost 
constant surface brightness with sharp edges. The above 
properties as well as the size and position of the bubbles 
are rather robust with respect to the details of foreground 
subtraction making it unlikely that the bubbles are an arte- 
fact of the foreground subtraction. 

Both the position at galactic longitude I = 0° and its sym- 
metry with respect to the galactic plane hint at the galactic 
centre (GC) as the origin of the bubbles. While similar 
structures have been observed in radio galaxies the detec- 
tion of the Fermi bubbles is puzzling given that there is 
no evidence for present activity of the massive black hole 
at the GC. Understanding this would provide an excellent 
probe of this region which is otherwise obscured by the 
galactic disk. The bubbles may play an important role in 
the dynamics of our galaxy and constitute a source of cos- 
mic rays. They dominate the high latitude 7-ray emission 
at (and possibly contribute close to) the GC, so constitute 
an important background for indirect dark matter searches. 
It is therefore important to understand and model the origin 
of the non-thermal emission from the bubbles. 



The observed high-energy gamma-rays can in principle be 
of hadronic or leptonic origin, i.e. tt° decay or inverse 
Compton (IC) scattering of the ambient radiation fields 
(CMB, far IR and optical/UV) by high energy electrons. A 
model [4] for the hadronic origin invokes an increased star 
formation rate close to the GC in combination with a strong 
convective wind up to kiloparsec distances from the galac- 
tic plane. The neutral pions are produced from spallation 
on the ambient gas in the bubbles of an E~ 2 spectrum of 
protons and nuclei accelerated in supernova remnants close 
to the GC; the kinematics of the ir° production can explain 
the spectral shoulder around 1 GeV. Leptonic models, on 
the other hand, have to explain how despite the rapid cool- 
ing of electrons the bubbles volume can be filled with a 
consistently hard spectrum. It has been argued that disrup- 
tion of stars close by the central massive black hole can heat 
the ambient gas and produce shocks. Estimates for the rate 
of this process predict hundreds of concentric shock fronts 
filling the bubbles. Electrons are repeatedly accelerated by 
diffusive shock acceleration, thereby explaining the hard 
gamma-ray spectrum. The spectral cut-off/roll-over would 
then be due to competition of acceleration and cooling by 
IC (and possibly synchrotron) losses. 

Data from the ROSAT x-ray satellite [5| however only 
show evidence for a limb brightened structure coinciding 
with the bubble edges, possibly from a shock front. The 
non-observation of x-rays from the bubble interior, on the 
other hand, points at a relatively thin, hot plasma. With 
estimates for a gas density of n ~ 10~ 2 cm~ 3 and a 
temperature of T ~ 2 keV, the total energy in hot gas 
erg (l]. Furthermore, assuming velocities 
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typical for shock fronts in the interstellar medium gives 
— 10 7 ([//1000kms~ 1 )yr for the age of the bubbles at a 
projected distance of 10 kpc. 

Other possible scenarios for the generation of the bubbles 
than those above include jets emanating from the central 
black hole. While we choose to remain agnostic about the 
origin of the bubbles itself we note that a shock might have 
been produced by such a jet active for a few million years. 
It has recently been shown |6 1 that a light but overpressured 
jet powered by ~ 10 % of the Eddington luminosity leads 
to a shock coincident with the bubble edge and in agree- 
ment with the overall bubble shape. In the following, we 
will explain the non-thermal emission from the bubbles by 
2nd-order Fermi acceleration of electrons and IC scattering 
of these electron on ambient radiation fields. As the elec- 
trons are constantly accelerated in the whole bubble this 
can lead to the hard gamma-ray spectrum. 
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Figure 1 : Relevant timescales as a function of energy. The 
acceleration and escape time depend on the distance £ = 
x/L from the shock and are shown close to the shock (£ = 
0.1) and in the bubble interior (£ = 1). 



2 Second order Fermi acceleration 

In particular, we start from the evidence for a shock front 
from ROSAT. At the outer shock Rayleigh-Taylor and 
Kelvin-Helmholtz instabilities will generate plasma turbu- 
lence that is then being convected into the bubble interior 
by the downstream bulk flow. The turbulence will cas- 
cade from the injection scale L to smaller scales and will 
finally be dissipated at a scale Id, i.e. once the eddy ve- 
locity reaches the Alfven velocity, w c dd('d) ~ va- The 
usual Rankine-Hugoniot conditions allow to compute 
the spatial variation of the eddy velocity at the injection 
scale, u, and the magnetosonic phase velocity, v-p, with dis- 
tance x = from the shock: 
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where U is the shock velocity, va the Alfven velocity 
(which we assume to be constant and equal to the speed 
of sound v St o at the shock) and a = 3 — I6v 2 /U 2 . 

We consider the stochastic acceleration by large-scale, fast 
mode turbulence [8 1. Second order Fermi acceleration pro- 
cesses like this have been proven successful in explaining 
the non-thermal spectra of high-energy electrons in a vari- 
ety of astrophysical environments |]9] [TOJ [TT] |T2j [T3] |7 1 and 
might be responsible for the acceleration of ultra-high en- 
ergy cosmic rays lfl4l H31 [l6l . The spectrum is governed 
by the Fokker-Planck equation, 
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where n(p, i) dp is the density of electrons with momentum 
in [p, p + dp] . The second, third and fourth term describe 
diffusion and systematic gains in momentum, escape due to 
spatial diffusion and energy losses by synchrotron radiation 



and IC scattering (with cooling time t coo \ ~ p/ (dp/ dt)), 
respectively. The diffusion coefficient in momentum for 
scattering by fast magnetosonic waves is JS) 
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which translates into the timescale for acceleration, i acc ~ 
p 2 /D pp . Both i acc and t csc (and therefore the resulting 
spectrum) depend on three parameters that cannot be in- 
ferred directly from observations. The scale of turbu- 
lence injection L is necessarily smaller than the size of 
the bubbles and MHD simulations show generation of tur- 
bulence on kiloparsec scales. Here, we assume L = 
2 kpc. The shock velocity can in principle be determined 
from the displacement of the shock; the shock needs ~ 
50 (<7/10 8 cms -1 ) yr to move a distance corresponding to 
the 1" resolution of the Chandra X-ray observatory. Here, 
we fix U — 2.6 x 10 s cms -1 , a value consistent with MHD 
simulations |6|. Finally the Alfven velocity is given by the 
square root of the ratio of magnetic field energy density to 
thermal plasma energy density: /3 a = «a/c = ^JUb/U p . 
Hence /3a > 2.8 x 10 -4 for an estimated upper limit on 
the thermal gas density n < 10 -2 cm -3 |T) and a mag- 
netic field B — 4 /xG. Such a field strength in the halo is 
suggested by radio observations of edge-on spiral galaxies 
as NGC 891 flT7). Here we adopt /3 A = 5 x 10 -4 . The 
resulting timescales are compared in Fig.Q] 

With these adopted parameters the dissipation length 
l d >8x 10 19 (L/kpc)(?7/10 8 cms- 1 )- 3 (^ A /10 -3 ) 3 cm 
is always larger than the gyro-radius of relativistic elec- 
trons - 7.5 x 10 n (B/4/iG) -1 (i;/GeV)cm. Hence, 
the spatial diffusion coefficient and also the escape and 
acceleration time are effectively energy-independent. 
Furthermore, with the parameters as above we recover a 
hierarchy of timescales, t acc ,t csc <C £iif e which justifies 



the use of the steady state solution [18], 
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The spectral index, — a = 1/2 — y9/4 + t acc /t csc , is de- 
termined by the ratio of acceleration and escape times alone 
and asymptotically approaches — 1 as t acc /t csc — > 0. One 
could argue that for low energies the cooling time becomes 
larger than the dynamical time scale t\n c and that therefore 
the use of the steady state solution is not strictly justified. 
However, it has been shown [19| for 2nd-order Fermi ac- 
celeration that irrespective of the cooling rate the spectrum 
always attains the steady state spectrum in a few times t a cc 
and that the steady state solution can therefore be applied 
as long as i acc <C tuf e . The change in timescales i acc and 
t csc with distance from the shock front £ (through u(£) and 
Wp(0) * s therefore adiabatic such that the electron spec- 
trum relaxes quickly to its steady state value. 

The emissivity in gamma-rays is calculated in the most 
general form ll20l using a recent model of interstellar ra- 
diation fields ED and depends on the distance from the 
shock. We calculate the flux of photons by integrating 
along the line of sight through the bubble. Since the overall 
normalisation of the electron spectrum depends on the mi- 
crophysics of injection, we fix this by demanding that our 
model matches the observed gamma-ray flux. 
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Figure 3: The overall spectrum E 2 J 1 in gamma-rays. The 
data are shown as obtained with two different IC tem- 
plates [ 1 1. The fits from a hadronic [4| and a leptonic [22] 
model are shown by the dashed and dotted lines, respec- 
tively. The gamma-ray flux from our model is shown by 
the solid line and the dot-dashed lines show the contribu- 
tions from IC scattering on the CMB, FIR and optical/UV 
(from left to right). 



3 Results 

Figure [2] shows the electron spectrum for different dis- 
tances from the shock. The spectrum is hardest close to the 
shock and becomes gradually softer towards the bubble in- 
terior. Furthermore, the spectral pile-up and cut-off (deter- 
mined by the competition between acceleration and energy 
losses) move to lower energies. Integration of the (position- 
dependent) spectrum over both bubbles shows that the total 
energy in electrons above 100 McV is ~ 10 51 erg. This 
is a rather moderate energy demand, in particular in com- 
parison to the hadronic model |4| which requires up to five 
orders of magnitude more energy in high energy protons. 

The overall spectrum of gamma-rays from the bubbles is 
shown in Fig. [3] and compared to the data and predictions 
from other models. Our model not only reproduces the 
E~ 2 spectrum but also both the spectral shoulder around 
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Figure 2: The electron spectrum E 2 n e - at different dis- 
tances £ = x/L from the shock front. 



1 MeV and the roll-over/cut-off at ~ 200 GeV. We note 
that the other two models presented would need to invoke 
a somewhat unmotivated break in the proton/electron spec- 
trum to produce both features whereas in our model they 
arise naturally due to the very hard electron spectrum and 
the cut-off due to cooling. 

In Fig. [4] we compare the data with the intensity as a func- 
tion of distance from the bubble edge predicted by our 
model and obtained in the same fashion as in Ref. Q), i.e. 
averaging over great circles intersecting the bubble centre. 
At both energies for which data is available (2 and 10 GeV) 
our model nicely reproduces the constant profile inside the 
bubbles and their sharp edges, i.e. the jump in intensity 
within a few degrees around the bubble edge. We have also 
computed the profile at 500 GeV which is much more limb- 
brightened - a robust prediction of our model. We note 
that the profile expected from a constant volume emissivity, 
as predicted by a hadronic [4| and a leptonic ll22l model, 
would be much softer at the edges and does not reproduce 
the data. 

Both the constant surface brightness with sharp edges at 
low energies and the limb-brightening at the highest ener- 
gies are consequences of the position dependent electron 
spectrum. As shown already in Fig. [2] the high-energy 
electrons are present only close to the shock front whereas 
low-energy electrons are distributed over the whole bubble 
volume. The IC emissivity at a given distance from the 
shock is in fact a convolution over a certain range in elec- 
tron energies: For GeV gamma-rays the emissivity profile 
is flat in the bubble interior and peaks at the bubble edge, 
which in projection leads to a flat surface brightness with 
sharp edges. At the highest energies however, only elec- 
trons of hundreds of GeV which are located close to the 
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Figure 4: The intensity E 2 J 1 in gamma-rays is shown 
as a function of the distance from the bubble edge at 
2 GeV (solid line), 10 GeV (dashed line) and 500 GeV 
(dot-dashed line), together with the data [1] from the av- 
eraged 1 — 2 and 2 — 5 GeV and the averaged 5 — 10 and 
10 — 20 GeV maps. We also show the profile expected from 
a constant volume emissivity (dotted line) which clearly 
does not reproduce the observed profile. 



shock can contribute which leads to a limb-brightening in 
projection. 

While the "WMAP haze" |23| has not been observed in 
polarised emission [24 1 and may be just an artefact of 
the template subtraction |25|, it has been proposed as a 
physical counterpart of the Fermi bubbles However 
as seen in Fig. [5] in our model the expected synchrotron 
flux in the middle of the bubble is of the required am- 
plitude only if the magnetic field is as strong as 15 /iG. 
For a 4 /iG field the synchrotron flux is significantly lower, 
1.6 x 10- 21 (^/GHz)-°- 2 ergcm- 2 s- 1 sr- 1 . 
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Figure 5: Radio and microwave flux from synchrotron 
emission of electrons in the bubble for a magnetic field of 
4 /iG (solid line) and 15 /iG (dashed line). The data point 
shows the amplitude of the "WMAP haze" |26| together 
with a range of spectral indices and the inverse triangles 
are upper limits obtained from the 408 MHz all-sky sur- 
vey C2) and the WMAP 23 GHz, 33 GHz bands It28l . 



The hadronic model predicts a detectable flux of neutrinos 
for the proposed Mediterranean km 3 neutrino telescope [4]. 
We stress however that the observed bubble profile already 
disfavours this model (as well as the leptonic DSA model) 
and instead favours our model with 2nd-order Fermi accel- 
eration of electrons. 
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